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Abstract A future e+e^ collider, such as the ILC or CLIC, 
would allow the Higgs sector to be probed with a precision 
significantly beyond that achievable at the High-Luminosity 
LHC. A central part of the Higgs programme at an e+e^ col¬ 
lider is the model-independent determination of the absolute 
Higgs couplings to fermions and to gauge bosons. Here the 
measurement of the e+e^ -> HZ Higgsstrahlung cross sec¬ 
tion, using the recoil mass technique, sets the absolute scale 
for all Higgs coupling measurements. Previous studies have 
considered (7(e+e^ —> HZ) with Z —where £ = e, p. 
In this paper it is shown for the first time that a near model- 
independent recoil mass technique can be extended to the 
hadronic decays of the Z boson. Because the branching ra¬ 
tio for Z —> qq is approximately ten times greater than for 
Z —> £+£^, this method is statistically more powerful than 
using the leptonic decays. For an integrated luminosity of 
500fb^^ at a centre-of-mass energy of y/s = 350GeV at 
CLIC, (J(e +e —HZ) can be measured to ±1.8% using the 
hadronic recoil mass technique. A similar precision is found 
for the ILC operating at = 350GeV. The centre-of-mass 
dependence of this measurement technique is discussed, ar¬ 
guing for the initial operation of a future linear collider at 
just above the top-pair production threshold. 


1 Introduction 

A future e+e^ collider, such as the Compact Linear Collider 
(CLIC) [1] or the International Linear Collider (ILC) [2], 
would be complementary to the Large Hadron Collider (LHC) 
and the High-Luminosity LHC (HL-LHC), providing tests 
of beyond the Standard Model (SM) physics scenarios through 
a broad programme of highly precise measurements. A cen¬ 
tral part of the linear collider physics programme is the pre¬ 
cise study of the properties of the Higgs boson. The LHC 


and HL-LHC will provide an impressive range of Higgs 
physics measurements, establishing the general properties 
of the Higgs boson, such as its mass and spin. The LHC will 
also provide measurements of the product of the Higgs pro¬ 
duction rate and Higgs decay branching fractions into dif¬ 
ferent final states. Current estimates suggest that ratios of 
couplings can be measured to 2 % - 7 % (depending on the 
final state) with 3000 fb^^ of data [3]. A number of recent 
studies, see for example [3, 4], have indicated that modifica¬ 
tions of the Higgs couplings due to beyond the SM (BSM) 
physics are almost always less than 10 % and can be as small 
as 1 % - 2 % in a number of models. 

An e+e^ collider would be a unique facility for precision 
Higgs physics [5, 6], providing measurements of the Higgs 
boson branching ratios that may be an order of magnitude 
more precise than those achievable at the HL-LHC. Such 
measurements may be necessary to reveal BSM effects in 
the Higgs sector. Moreover, an e+e^ collider provides the 
opportunity to make a number of unique measurements in¬ 
cluding; i) absolute measurements of Higgs couplings, rather 
than ratios; ii) a precise measurement of possible decays to 
invisible (long-lived neutral) final states; and iii) a < 10% 
measurement of the total Higgs decay width, Th. In addi¬ 
tion, an e+e^ collider operating at 1 TeV or above, for ex¬ 
ample CLIC or an upgraded ILC, would have sensitivity to 
the Yukawa coupling of the top quark to the Higgs boson 
and the Higgs self-coupling parameter A, thus providing a 
direct probe of the Higgs potential. 

This paper presents the first detailed study of the poten¬ 
tial of making a model-independent measurement of guzz 
from the recoil mass distribution in e+e^ —HZ with Z ^ 
qq, denoted as HZ(Z —> qq). The studies were initially per¬ 
formed in the context of the CLIC accelerator operating at 
-v/s = 350GeV. The studies were repeated for the ILC oper- 
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ating at the same energy and for CLIC at y/s = 250 GeV and 
v^ = 420GeV. 


1.1 Higgs production in e+e collisions 


In e+e^ collisions at ^/s = 250 - 500 GeV, the two main 
Higgs production mechanisms are the Higgsstrahlung pro¬ 
cess, e+e^ —> HZ, and the WW-fusion process, e+e^ —> 
HVgVe, shown in Fig. 1. For oth 125 GeV, the cross section 
for the i-channel Higgsstrahlung process is maximal close to 
= 250 GeV, whereas the cross section for the f-channel 
WW-fusion process increases with centre-of-mass energy, 
as indicated in Tab. 1. 


e+ Z e+ Ve 




Fig. 1 Feynman diagrams for: (left) the Higgsstrahlung process 
e+e“ —> HZ, which dominates at yT = 250 GeV; and (right) the WW- 
fusion process e+e“ —> Vj.VeH, which dominates at yT > 500 GeV. 


The total HZ cross section is proportional to the square of 
the coupling between the Higgs and Z bosons, gnzz, 

(7(e+e- ^ HZ) oc ^2 , 

and the cross sections for the exclusive final-state decays 
H —XX can be expressed as 


(7(e+e^ HZ) X BR(H ^ XX) oc 
(7(e+e^ ^ Hv^Ve) x BR(H ^ XX) oc 


^Hzz X 8 hxx 


8rww X giixx 
Ui 


Once gHZZ has been determined in a model-independent 
manner, the ratio of the Higgsstrahlung and WW-fusion cross 
sections for the same exclusive Higgs boson final state yields 
gHWW- Subsequently, the measurement of (j(e+e^ -> Hv^Ve) 

X BR(H —^ WW*), which depends on gnww/^’ provides 
a determination of Th. At this point all measurements of ex¬ 
clusive Higgs decays provide absolute and model-independent 
determinations of the relevant coupling(s). The determina¬ 
tion of gHZZ from the recoil mass distribution in e+e ^HZ 
lies at the heart of this scientific programme. 


1.2 The leptonic recoil mass measurement 

The Higgsstrahlung process provides the opportunity to study 
the couplings of the Higgs boson in a model-independent 
manner. This is unique to an electron-positron collider, where 
the clean experimental environment and the relatively low 
SM cross sections for background processes allow e+e^ —>^ 
HZ events to be selected based solely on the measurement 
of the four-momentum of the Z boson, regardless of how 
the Higgs boson decays. The clearest topologies occur for 
Z —> e+e^ and Z —> decays, which can be identified by 

first requiring that the measured di-lepton invariant mass ma 
is consistent with mz - The four-momentum of the system re¬ 
coiling against the Z boson is obtained from Free = — £« 

and prec = —Ptt- In e+e^ —HZ events, the invariant mass 
of this recoiling system, mrec, will peak at mn. Fig. 2 shows 
the simulated recoil mass distribution in the ILD [8] detec¬ 
tor concept for 250fb^' of ILC data at y^s = 250GeV with 
beam polarisation P(e^,e+) = (—0.8, -1-0.3); by combining 
both Z —and Z —>■ e+e^ decays, (7(HZ) can be mea¬ 
sured to 2.6 %, leading to a determination of gnzz with a 
precision of 1.3 % [5, 8, 9]. 




rriH /GeV 

Fig. 2 The recoil mass distribution for the Higgsstrahlung process for 
a) Z —> |rX|x“ andb) Z —> e+e“ at yT = 250GeV. The results are shown 
for /’(e^,e+) = (—80%,-|-30%) beam polarisation. Taken from [8]. 
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e+e polarisation = 

250 GeV 

350 GeV 

500 GeV 

C7(e+e —>■ HZ) 

unpolarised 

211fb 

134 fb 

65 fb 

^(e+e^ —> HVgVe) 

unpolarised 

21fb 

34 fb 

72 fb 

^(e+e- HZ) 

(-0.8,+0.3) 

318fb 

198 fb 

96 fb 

^(e+e^ —>■ HVgVe) 

(-0.8,+0.3) 

37 fb 

73 fb 

163 fb 


Table 1 The leading-order Higgs cross sections for the Higgsstrahlung and WW-fusion processes for me = 125GeV at three centre-of-mass 
energies. The cross sections are calculated [7] including initial-state radiation and are shown for unpolarised electron/positron beams and assuming 
the baseline ILC polarisation of P(e“,e+) = (—0.8, -fO.3). 


1.3 Recoil mass measurement at different centre-of-mass 
energies 

The narrowness of the recoil mass peak is an important fac¬ 
tor in determining the precision to which gnzz can be mea¬ 
sured. The recoil mass can be expressed as 

™rec = EuY ~Pee = — Pa 

= S — 2y/s {El\ +E(2) + tWa : 

where is the centre-of-mass energy and E^ and Ee 2 are 
the energies of the two leptons. Since ma will peak narrowly 
around mz, it can be seen that the width of the recoil mass 
peak scales with both and the lepton energy (or momen¬ 
tum) resolution. For high-momenta muons, where multiple 
scattering in the tracking chambers is relatively unimpor¬ 
tant, the fractional momentum resolution Gp/p will scale 
approximately as the transverse momentum pj, thus 
will scale quadratically as p - pj. Consequently, in the range 
^/s = 250 — 500 GeV, where the energy of the fermions from 
the Z decay approximately scales as ^/s, the width of the 
recoil mass distribution increases significantly with increas¬ 
ing centre-of-mass energy. For this reason, the leptonic re¬ 
coil mass analysis leads to a higher precision on gazz for 
~ 250 GeV [5], where the ff (HZ) is largest and the re¬ 
constructed recoil mass peak is relatively narrow. This has 
been one of the strongest arguments for the initial operation 
of the ILC at a relatively low centre-of-mass energy. This ar¬ 
gument does not apply to the recoil mass measurement with 
hadronic Z —>^ qq decays since the recoil mass resolution 
depends less strongly on ^/s than for leptonic final states be¬ 
cause the jet energy resolution for the linear collider detector 
concepts scales linearly with energy, Og ~ O.OSF [10]. Al¬ 
though the hadronic recoil mass measurement has been con¬ 
sidered previously [11], this paper presents the first detailed 
study of its potential. 


2 Monte Carlo Samples, Detector Simulation and Event 
Reconstruction 

The CLIC results presented in this paper are based on de¬ 
tailed Monte Carlo (MC) simulation using: a full set of SM 


background processes; a detailed Geant 4 [12, 13] simu¬ 
lation of the CLICJLD detector concept [14]; and a full 
reconstruction of the simulated events. 


2.1 Monte Carlo event generation 

The simulated SM event samples were generated using the 
WHIZARD 1.95 [7] program. The expected energy spec¬ 
tra for the CLIC beams, including the effects from beam- 
strahlung and the intrinsic machine energy spread, were used 
for the initial-state electrons and positrons. The process of 
fragmentation and hadronisation of final-state quarks and 
gluons was simulated using PYTHIA 6.4 [15] with a pa¬ 
rameter set [16] that was tuned to OPAL e+e^ data recorded 
at LEP. The decays of T leptons were simulated using the 
TAUOLA package [17]. The mass of the Higgs boson was 
taken to be oth = 126 GeV and the decays of the Higgs bo¬ 
son were simulated using PYTHIA with the branching frac¬ 
tions of [18]. A dedicated sample of e^e —^ HZ events with 
Higgs decays to “invisible” long-lived neutral particles was 
produced by artificially setting the Higgs boson lifetime to 
infinity. Because of the 0.5 ns bunch spacing in the CLIC 
beams, the pile-up of beam-induced backgrounds from the 
YY —^ hadrons process was included in the simulated event 
samples to ensure its effect on the event reconstruction was 
accounted for. 


2.2 CLIC detector simulation and event reconstruction 

The GEANT4-based Mokka [19] program was used to sim¬ 
ulate the detector response of the CLICJLD detector con¬ 
cept [14]. The QGSP_BERT physics list was used to model 
the hadronic interactions of particles in the detectors. The 
hit digitisation and the event reconstruction were performed 
using the Marlin [20] software packages. Particle flow re¬ 
construction was performed using PandoraPFA [10, 21]. An 
algorithm, using the individual reconstructed particles, was 
used to identify and remove approximately 90 % of the out- 
of-time background due to pile-up from yy hadrons; here 
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the Loose particle flow object selection, described in [14], 
was used. 

Jet finding was performed using the FastJet [22] package. 
Because of the presence of pile-up from yy —hadrons, the 
ee_kt (Durham) algorithm employed at LEP is not effective 
as it clusters particles from pile-up into the reconstructed 
jets. Instead, the hadron-collider inspired k, algorithm, with 
the distance parameter R based on A rj and A(j), was used 
with R = n/2. This algorithm allows particles to be clustered 
into “beam jets”, aligned with the beam axis, in addition to 
jets seeded by high-momentum particles. Background from 
the pile-up of yy —^ hadrons can, to a large degree, be re¬ 
moved by ignoring particles in the “beam jets”, largely mit¬ 
igating the impact of the beam background. 

The hadronic recoil mass study, presented in this paper, cov¬ 
ers a wide range of HZ(Z —> qq) final-state topologies rang¬ 
ing from two jets where Higgs decays to long-lived neutral 
particles, e+e^ —)• HZ (invis.)(qq), to six-jet toplogies 
from, for example, e+e^ —^ HZ —(WW* —qqqq)(qq). 
For this reason, each reconstructed event is clustered into 
two-, three-, four-, five- and six-jet topologies, with “y-cut” 
variables used to indicate the underlying physical topology. 
For example, if an event is forced into a three-jet topology, 
y 34 is the k, value at which the event would be reconstructed 
as four jets and y 23 is the k, value at which the event would 
be reconstructed as two jets. 

2.3 ILC detector simulation and event reconstruction 

The event generation and reconstruction for the ILC studies, 
presented in section 4, follows closely that described above. 
The main differences are: i) the ILC beam spectrum, where 
the effects of beamstrahlung are less pronounced, was used, 
ii) the detector simulation used the ILD detector concept for 
the ILC, rather than the CLIC_ILD model adapted for CLIC; 
and iii) the much longer bunch spacing at the ILC means 
that only in-time background from yy —> hadrons needs to 
be included. 

3 Hadronic Recoil Mass Measurement at CLIC 

In the process e+e^ -> HZ it is possible to cleanly identify 
Z —>^ e+e^ and Z —> decays regardless of the H decay 

mode. Consequently, the selection efficiency is almost inde¬ 
pendent of the nature of the H decay. For Z —> qq decays, the 
selection efficiency will depend more strongly on the Higgs 
decay mode. For example, in (H —!• bb)(Z —> qq) events, 
the reconstruction of the Z boson is complicated by mis- 
associations of particles to jets and by the three-fold ambigu¬ 
ity in associating four jets to the Z and H. These ambiguities 


will increase with the number of jets in the final state. For 
this reason, it is much more difficult to construct an event se¬ 
lection, based only on the reconstructed candidate Z —> qq 
decay, with a selection efficiency that is independent of the 
Higgs decay mode. Nevertheless it is possible to minimise 
this dependence. The strategy adopted here is to: i) sepa¬ 
rate all simulated events into candidates for Higgs decays to 
“invisible” long-lived neutral particles and decays to visible 
final states; ii) identify the di-jet system that is the best can¬ 
didate for the Z ^ qq decay; iii) reject events consistent with 
a number of clear background topologies using the informa¬ 
tion from the whole event; iv) identify HZ(Z —qq) events 
solely based on the properties from the candidate Z —qq 
decay, first for the candidate visible Higgs decays and then 
for the candidate invisible Higgs decays; and v) combine the 
results into a single measurement of ( 7 (HZ). 


3.1 Separation into candidate visible and invisible Higgs 
decay samples 

Hadronic events are selected by forcing each event into a 
two-jet topology and requiring at least three charged parti¬ 
cles in each jet. The surviving events are then divided into 
candidates for either visible H decays or invisible H decays, 
in both cases produced in association with a Z —> qq. Events 
are categorised as potential invisible H decays on the ba¬ 
sis of the y-cut values in the k, jet-finding algorithm. For 
invisible H decays, only the Z —> qq system is visible in 
the detector, typically resulting in a two-jet topology (with 
the possibility that QCD radiation can increase the number 
of reconstructed jets). Consequently, invisible H decays will 
have small values of y 23 and ^ 34 , the variables respectively 
representing the k, value at which an event transitions from 
two to three jets and from three to four jets, as indicated 
in Fig. 3. Events are categorised as candidate invisible H 
decays if —logio(y 23 ) > 2.0 and —logio(y 34 ) > 3.0. Due 
to gluon radiation in the parton shower, only 74 % of the 
simulated HZ(Z —?> qq) events with invisible H decays are 
placed in this two-jet topology candidate invisible H decay 
sample. To improve the efficiency for correctly categoris¬ 
ing SM Higgs decays with low-energy leptons, for exam¬ 
ple H WW* — xvTV, events with — logjo(y 23 ) < 2.5 and 
— Iogjo(y 34 ) <3.5 are forced into three jets and are excluded 
from the invisible Higgs decay sample if the lowest-energy 
jet has fewer than four reconstructed tracks or contains an 
identified e*/p* with energy E > 5GeV. Only 2.2% of 
simulated HZ events with SM Higgs decays end up in the 
candidate invisible Higgs sample. 
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Fig. 3 The distributions of — log[o(>' 23 ) and — logjo(>’ 34 ) for simulated HZ(Z —t qq) events for visible and invisible Higgs decays at = 
350GeV. The distributions are normalised to an integrated luminosity of 500 fb“*. The distribution for the invisible H decays assumes a 100% 
branching fraction into invisible decay modes. The vertical lines with arrows indicate the cut values used in this analysis. 


3.2 Recoil mass reconstruction 

For each candidate HZ(Z qq) event, the recoil mass is 
calculated from = {^/s — Eqq )^ ~Pqq’ where Fqq and 
Pqq are the summed energy and momentum of the di-jet sys¬ 
tem from the identihed candidate Z —> qq decay. In the case 
of the candidate invisible Higgs decay sample, the two jets 
are assumed to be from Z —qq. The resulting recoil mass 
distribution for candidate invisible Higgs decays, which is 
strongly peaked around tWrec is shown in Fig. 4a. In 

the case of the candidate visible Higgs decay sample, the 
situation is more complicated as this sample encompasses 
many different HZ event topologies. For example, H -3- bb 
decays will result in a four-quark HZ hnal state, usually 
yielding four jets, whereas, H —WW* —qqfv and H —^ 
WW* —7> qqqq decays will respectively usually yield hve- 
and six-jet hnal states. In all cases gluon radiation in the 
parton shower can increase the reconstructed jet multiplic¬ 
ity relative to the tree-level expectation. 

In order to achieve the desired (near) model independence 
of the analysis, it is necessary to have a similar quality of 
recoil mass reconstruction for all Higgs boson visible decay 
modes. This hinges on the correct identihcation and recon¬ 
struction of the Z —> qq di-jet system. The hrst stage is to 
force events in the candidate visible Higgs decay sample into 
a four-jet topology. From the three possible di-jet combina¬ 
tions, the di-jet system with invariant mass mqq closest to mz 
is identihed as the candidate Z —> qq decay and its energy 
and momentum are used to calculate the recoil mass mrec- In 
selecting the candidate Z decay, only jets containing more 
than three charged particles are considered. To improve the 
reconstruction of higher-jet-multiplicity hnal states, such as 


WW* —qqqq, each event is also forced into hve jets 
and the di-jet system with mass closest to mz is again iden¬ 
tihed as the candidate Z —qq decay. The hve-jet topology 
is only used if — logio(y 45 ) >3.5 and both Wqq and miec are 
respectively closer to mz and mu than the corresponding val¬ 
ues from the four-jet reconstruction. Even in the genuine six- 
parton topology HZ —^ (WW*)qq —^ (qqqq)qq only 13% 
of events are reconstructed as hve jets, for the remainder, the 
four-jet reconstruction is preferred. However, provided the 
jets from the Z —> qq decay are correctly identihed, there is 
no need to correctly reconstruct the recoiling system as only 
the properties of the Z —> qq decay are used in the subse¬ 
quent analysis. For this reason, allowing the possibility of 
reconstructing events as six jets was found not to improve 
the overall recoil mass reconstruction. Fig. 4b shows the re¬ 
sulting recoil mass distribution for simulated HZ events with 
H bb, H —>^ WW* —?► qqqq and H —> Despite the 

very different hnal states, similar recoil mass distributions 
are obtained. 


3.3 Preselection 

After dividing all events into either candidates for visible 
or invisible Higgs decays and having identihed the two jets 
forming the candidate Z —qq system, preselection cuts are 
applied to reduce backgrounds from larger cross section SM 
processes such as e+e^ —> qq and e+e^ —> qqqq. Cuts are 
based on the invariant mass of the Z —qq candidate, mqq, 
and corresponding recoil mass, mrec- In addition, the invari¬ 
ant mass of all the visible particles not originating from the 
candidate Z —^ qq decay, m^, is calculated. It is important 
to note that m% is only used to reject specihc background 
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Fig. 4 a) The reconstructed hadronic recoil mass distribution for the candidate HZ events with Z —> qq and H —> invis. b) The reconstructed 
hadronic recoil mass distributions for candidate HZ(Z —> qq) and either H —> bb, H —> WW* —> qqqq or H —> In each case the distributions 

are normalised to unit area. An underflow (not shown) contains the small fraction of events where no good Z —> qq candidate is identified. 


topologies in the preselection and is not used in the main se¬ 
lection; in HZ events will depend strongly on the Higgs 
decay mode. The preselection cuts (most of which are com¬ 
mon to the visible and invisible Higgs selections) are: 

- the event must be broadly consistent with being HZ, 
70 GeV < mqq < 110 GeV and 80 GeV < Wrec < 200 GeV. 

- background from e+e^ —^ qq is suppressed by removing 
events with net transverse momentum pj < 20 GeV and 
— Iogio(y 34 ) > 2.5, indicating a final-state system con¬ 
sisting of fewer than four primary particles. 

- events in the invisible Higgs decay sample are rejected 
if |cos0jnis| > 0.7, where is the polar angle of the 
missing momentum vector, almost completely eliminat¬ 
ing the contribution from e'^e —^ qq with unobserved 
initial-state radiation (ISR). 

- events in the invisible Higgs decay sample are rejected 
if there is an isolated identified e*/p* with energy > 
lOGeV, suppressing background from WW —> qqfv. 

- background from e+e^ —> qq with unobserved ISR, in¬ 
cluding radiative return to the Z resonance, is suppressed 
by rejecting events with net transverse momentum px < 
20GeV and |cos0n,i,s| > 0.9. 

- the background from e+e^ ^ WW —> qqqq is suppressed 
by forcing events into four jets and selecting the di-jet 
pair with the mass mqq closest to mw- Events are re¬ 
jected if Pi < 20GeV and 65GeV < mqq < lOOGeV 
and 65 GeV < mV < 100 GeV, where mV is the mea¬ 
sured invariant mass of the second di-jet pair. 

- the background from e+e^ —ZZ —qqqq is suppressed 
in a similar manner. Events are forced into four jets and 
the di-jet pair with the mqq closest to mz is identified. 
Events are rejected if pj < 20GeV and 70GeV < mqq < 


105 GeV and 70 GeV < mV < 105 GeV, where mV is 
the measured mass of the second di-jet pair. 

The effects of the preselection cuts are summarised in Tab. 2. 
The events passing the preselection cuts are put forward as 
candidate HZ(Z —> qq) events with either: i) visible H de¬ 
cay products; or ii) invisible H decay products, depending 
on whether the event was consistent with a two-jet topology 
or not. The first two cuts listed above result in the largest 
loss of signal efficiency for the visible Higgs decay selec¬ 
tion. The qq£v background in the visible Higgs decay pre¬ 
selection could have been significantly reduced by reject¬ 
ing events with visible high-energy isolated leptons, but this 
would have introduced a bias against H decays with leptons 
in the final state. 


Final state 

ff/fb 

pvis. 

^presel 

pinvis. 

*^presel 

"plesel 

^invis 

presel 

qq 

25180 

0.4% 

— 

54570 

0 

qqfv 

5914 

11.2% 

0.9% 

326420 

26060 

qqqq 

5847 

3.8% 

- 

110520 

0 

qqU 

1704 

1.5% 

- 

13260 

60 

qqvv 

325 

0.6% 

14.8 % 

1050 

24180 

Hv,v, 

52 

2.5% 

5.6% 

640 

1430 

HZ; Z —> qq 

93 

42.0 % 

0.2% 

19630 

80 

H-> invis. (100%) 93 

0.6% 

48.6 % 

300 

22710 


Table 2 Summary of the effects of the preselection cuts for the the vis¬ 
ible and invisible recoil mass analyses. The efficiencies e include the 
effects of the preselection cuts and the division into the candidate visi¬ 
ble and invisible Higgs decay samples. The expected numbers of events 
passing the preselection cuts correspond to an integrated luminosity of 
500fb“* at CLIC, assuming unpolarised beams at yT = 350GeV. The 
numbers shown for the invisible Higgs decay modes correspond to a 
100 % branching ratio. 
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3.4 Selection of HZ-^qq with visible Higgs decays 


After preselection, the main backgrounds in the visible Higgs 
decay analysis arise from e+e^ —> qqfv and e+e^ —> qqqq, 
dominated by WW, single-W (Wev^) and ZZ processes. The 
event selection is based entirely on the reconstructed can¬ 
didate Z —> qq decay in the event. The properties of the 
remainder of the event (or the event as a whole) are not 
used as their inclusion would break the desired model in¬ 
dependence of the selection. For example, the background 
from e+e^ —qqfv can be significantly reduced by placing 
a lower bound on the total visible energy in the event, how¬ 
ever such a cut would bias the selection against Higgs decays 
with missing energy, such as H x+t^. 


The event selection uses a relative likelihood approach with 
discriminant variables based on the properties of candidate 
Z —^ qq decay. Two event categories are considered: a) the 
e+e^ —> HZ —> Hqq signal; and b) all non-Higgs background 
processes. The relative likelihood for an event being classi¬ 
fied as signal is defined as 


if = 


Tsignal 



where the individual absolute likelihood Lj for the event 
class j (signal or background) is formed from normalised 
probability distributions Pj{xi) of the discriminant variables 
Xi for that event class j: 


— ^preset ^ 11 ’ 

/ 


where CJpj.gj,g[ is the cross section after preselection for event 
class j. 


The discriminant variables used in the likelihood selection, 
all of which are based on the candidate Z —> qq decay, are: 
i) the two-dimensional distribution of mqq and ii) the 
polar angle of the Z candidate, |cosfe|; and iii) the mod¬ 
ulus of the angle of the jets from the Z —> qq decay in Z 
rest frame, relative to its laboratory frame direction of mo¬ 
tion, I cos 0q |. The two-dimensional distributions of mqq and 
mrec, are shown separately for the signal and background 
in Fig. 5. As expected, the HZ signal events peak around 
ffjqq « OTz and mrec ~ The anti-correlation between mi-ec 
and mqq is expected; when the reconstructed jet energies are 
higher than the true energies, the reconstructed value of mqq 
will be higher than mz and mrec will be lower than mH due 
to the —2y/sEz term in the expression for the recoil mass, 
m^gc = s — 2^/sEx + m^q . The broad peaked structure in the 
background distribution at lower values of mqq arises from 
qqfv events (which have been forced into a four- or five-jet 
topology). The use of the two-dimensional distribution of 
mqq versus mrec in the likelihood accounts for the associated 
correlations. 
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Fig. 5 The distribution of the reconstructed Z —> qq mass, mqq, versus 
the hadronic recoil mass, mrec, for HZ(Z —> qq) events (top) and for 
all background processes (bottom). In both cases the plots show all 
events passing the visible Higgs preselection for CLIC operating at 
= 350GeV. 




The two angular variables used in the likelihood selection 
are shown in Fig. 6. The discriminating power arises from 
the fact that the Higgs boson is a scalar particle, and the an¬ 
gular distributions in HZ production are different from those 
in the dominant backgrounds which mostly arise from the 
production of two vector particles. 

The resulting relative likelihood distribution is shown in Fig. 7 
Despite the fact that the signal-to-background ratio in the 
preselected event sample is approximately 1:25, the likeli¬ 
hood selection provides good separation. The statistical pre- 
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Fig. 6 a) The polar angle of the reconstructed Z candidates, | cos ^ |, for both signal and background events for CLIC operating at = 350GeV, 
and b) the modulus of the angle of the jets from the Z —> qq decay relative to the Z direction after boosting into its rest frame, | cos 0q |. The signal 
and background distributions are normalised to 500 fb“*, but the HZ(Z —> qq) signal has been scaled by a factor of 25 to improve its visibility. 


cision on the cross section for HZ production (where the 
Z decays hadronically and the H has SM branching frac¬ 
tions) is maximised with a likelihood cut of .if > 0.65. The 
resulting efficiencies and the expected numbers of selected 
events for an integrated luminosity of 500 fb^^ are shown 
in Tab. 3. The corresponding statistical uncertainty on the 
production cross section is ±1.9%. The precision can be 
improved by extracting the number of signal events by per¬ 
forming a maximum-likelihood ht to the shape of the sim¬ 
ulated likelihood distribution by varying the normalisations 
of the signal and background components, yielding a statis¬ 
tical error of 

^tTvis. = ±1.7% . 


Process 

cr/fb 

^resel 

c.vis. 

^^>0.65 

N^>Q.65 

qq 

25180 

0.4% 

0.07 % 

8525 

qqfv 

5914 

11.2% 

0.20 % 

5767 

qqqq 

5847 

3.8% 

0.49 % 

14142 

qqU 

1704 

1.5% 

0.22 % 

1961 

qqvv 

325 

0.6% 

0.04 % 

60 

Hv^Ve 

52 

2.5% 

0.23 % 

60 

HZ; Z —> qq 

93 

42.0% 

22.6 % 

10568 

H —> invis. (100%) 

[93] 

0.6% 

0.04 % 

20 


Table 3 Summary of the CLIC (H —> vis.)(Z —> qq) event selection at 
vT = 350 GeV, giving the cross sections, preselection efficiency, over¬ 
all selection efficiency for a likelihood cut of > 0.65 and the ex¬ 
pected numbers of events passing the event selection for an integrated 
luminosity of 500 fb“* assuming unpolarised electron and positron 
beams. The numbers shown for the invisible Higgs decay modes cor¬ 
respond to a 100 % branching ratio. 


c 

^ 30000 
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Fig. 7 The resulting likelihood distribution for the hadronic recoil 
mass analysis. The distributions correspond to 500 fb“* of CLIC oper¬ 
ation at yT = 350 GeV with unpolarised electron and positron beams. 
The optimal likelihood cut at .if = 0.65 is indicated by the arrow. 


-H(Z^qq) 

—qqlv 

■qqqq 

■qq 

■ other background 



3.5 Selection of HZ —5 qq with invisible Higgs decays 

The main backgrounds after preselection for the invisible 
Higgs decay selection arise from e+e^ —> qqfv and e+e^ —5 
qqvv, which are dominated respectively by single-W (WeVg) 
and ZZ processes. A relative likelihood selection is used 
to separate the (H —5 invis.)(Z —5 qq) signal from the non- 
Higgs background. The discriminant variables employed are 
the same as those used for the visible Higgs decay likeli- 
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Process 

(j/fb 

^presel 

c-invis. 

^^>0.60 

Vif>0.60 

qqiv 

5914 

0.9% 

0.03 % 

951 

qqvv 

325 

14.8 % 

1.83% 

2985 

Hv^Ve 

52 

5.6% 

0.37 % 

95 

HZ; Z —> qq 

93 

0.2% 

0.06 % 

31 

H—> invis. (100%) 

[93] 

48.6 % 

23.52% 

10983 


Table 4 Summary of the CLIC (H —> invis.)(Z —> qq) event selection 
at = 350 GeV, giving the raw cross sections, preselection efficiency, 
overall selection efficiency for a likelihood cut of ^ > 0.60 and the ex¬ 
pected numbers of events passing the event selection for an integrated 
luminosity of 500fb“* and unpolarised electron and positron beams. 
The numbers shown for the invisible Higgs decay modes correspond 
to a 100 % branching ratio. 


hood function, namely the two-dimensional distribution of 
(Wqq versus mjec, |cos0z| and |cos0q|. The most powerful 
of these variables is the recoil mass itself, shown in Fig. 8a, 
where the signal is plotted for the artificial case of BR(H 
invis.) = 100%. The resulting relative likelihood distribu¬ 
tion for signal and background is shown in Fig. 8b, where 
good separation between signal and background is achieved. 

In the limit where the H -9- invis. branching ratio is small 
(as expected), the expected uncertainty on the number of 
invisible Higgs decays selected by a particular likelihood 
cut is driven by the statistical fluctuations on the number 
of background events, v/B- In this limit, the correspond¬ 
ing uncertainty on the cross section for HZ production with 
H —^ invis. is given by 


An- ■ — 

^Omvis. — g <->^2 I 

where S is the number of signal events that would have been 
selected for the case of a 100 % branching fraction for H —> 
invis. This uncertainty is minimised for a relative likelihood 
cut of .if > 0.60, resulting in a ±0.58% statistical uncer¬ 
tainty on a Oinvi,s., relative to the SM cross section for e+e^ — 
HZ. The corresponding selection efficiencies are shown in 
Tab. 4, where the expected background from SM HZ pro¬ 
duction includes the H —ZZ* —^ vvvv component that has 
a SM branching fraction of 0.1 %. 

A more optimal approach to extracting the signal cross sec¬ 
tion is to fit the shape of the likelihood distribution of Fig. 8b, 
rather than simply imposing a single likelihood cut. In the 
limit that the invisible branching ratio is small, the result¬ 
ing Gaussian uncertainty on the HZ production cross section 
with invisible Higgs decays is 


A 


T.SM 

-'hz 


±0.56%, 


relative to the SM e+e —HZ cross section. For the SM 
Higgs, the corresponding expected 90% confidence level 


upper limit on the invisible Higgs branching ratio is 
BR(H ^ invis.) < 0.94% at 90% C.L. 


3.6 Model independence of the hadronic recoil mass 
measurement 


By combining the two analyses for HZ production where 
Z —> qq and the Higgs decays either to visible or invisible 
final states. 


ff(HZ) 


t^vis. ± t^invis. 

BR(Z ^ qq) ’ 


it is possible to determine the absolute e+e^ —> HZ cross 
section in a nearly model-independent manner. Since the 
fractional uncertainties (relative to the total cross section) on 
the visible and invisible cross sections are 1.7% and 0.6% 
respectively, the fractional uncertainty on the total cross sec¬ 
tion will be the quadrature sum of these two fractional un¬ 
certainties, namely 


A(7(HZ) = ±1.8%. 


Thus, the Higgsstrahlung cross section can be measured with 
a precision of better than 2% at y/s = 350 GeV using the 
hadronic recoil mass (with unpolarised beams). Such a mea¬ 
surement is competitive with that obtainable from the lep- 
tonic recoil mass measurement at = 250 GeV, where a 
precision of ±2.6%[5] is achievable with 250fb^' of data 
(assuming —80% and ±30% polarisation of the electron 
and positron beams). The strongest physics argument for 
operating a linear collider at = 250 GeV is the model- 

independent measurement of ff(HZ) that provides a deter¬ 
mination of gHZZ- If it can be argued that the hadronic recoil 
mass measurement is effectively independent of the nature 
of the Higgs boson decay modes (including possible exten¬ 
sions to the SM), then the arguments for operating an e+e^ 
linear collider at y/s ^ 250 GeV are greatly reduced; almost 
all other measurements of the properties of the Higgs boson 
are found to benefit from higher centre-of-mass energies [5]. 
In addition, operating at y/s ^ 350 GeV allows the study of 
Higgs production through the WW-fusion process and the 
pair production of top quarks. 

However, the hadronic recoil mass measurement of (j(HZ) 
can only be truly model independent if the overall (visible 
+ invisible) selection efficiency is independent of the Higgs 
decay mode. Tab. 5 summarises the combined selection ef¬ 
ficiency for e+e^ —!• HZ(Z —qq), broken down into the 
different Higgs decay modes. Also shown are the efficiency 
for H —> WW* decays broken down into the different W de¬ 
cay modes, covering a very wide range of event topologies. 
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Fig. 8 a) The reconstructed hadronic recoil mass distribution for events passing the preselection cuts in the clear two-jet topology, b) The invisible 
Higgs decay relative likelihood distribution for signal and background. In both distributions the event rates are normalised to a CLIC integrated 
luminosity of 500fb“* at •/s — 350GeV. The HZ signal is shown for the artificial case of BR(H —> invis.) = 100%. The optimal likelihood cut at 
Jif = 0.60 is indicated by the arrow. 


Decay mode 

cVlS. 

*^if>0.65 

cinvis. 

gVis. _|_ ginvis. 

H —> invis. 

<0.1% 

23.5 % 

23.5 % 

H-i-qq/gg 

22.6 % 

<0.1 % 

22.6 % 

H-> WW* 

22.1% 

0.1 % 

22.2 % 

H-l-ZZ* 

20.6 % 

1.1% 

21.7% 

H —>■ 

25.3 % 

0.2% 

25.5 % 

H —> yy 

25.7 % 

<0.1 % 

25.7 % 

H->Zy 

18.6% 

0.3% 

18.9% 

H —> WW* —)■ qqqq 

20.8 % 

<0.1 % 

20.8 % 

H -> WW* qqfv 

23.3 % 

<0.1 % 

23.3 % 

H —> WW* —> qqxv 

23.1% 

<0.1 % 

23.1 % 

H WW* -> fvfv 

26.5 % 

0.1% 

26.5 % 

H WW* -> ivxv 

21.1% 

0.5% 

21.6% 

H —>■ WW* —)■ xvxv 

16.3% 

2.3% 

18.7% 


Table 5 Summary of the efficiencies of the H(Z —> qq) analyses at 
= 350GeV, giving the overall selection efficiency for the visible 
analysis (.5f > 0.65) and the invisible Higgs analysis (.5? > 0.60). Here 
£ refers to either e or p. 

from four-jet final states (qqqq) to final states with two rela¬ 
tively soft particles, for example the visible tau decay prod¬ 
ucts from WW* xvrv. For all final-state topologies, the 
combined (visible + invisible) selection efficiency lies be¬ 
tween 19 % and 26 % compated to the mean selection effi¬ 
ciency of ^23%; a relative variation of ±15%. It should 
be noted that these numbers are only indicative, since the 
measured cross sections are extracted from fits to the like¬ 
lihood distributions, rather than from a selection imposing 
hard cuts. 

To assess the impact of the different sensitivities to the dif¬ 
ferent H decay topologies, the different Higgs decay modes 
in the HZ MC samples are reweighted to correspond to mod¬ 


Decay mode 

4(BR) 

O-vis. pinvis. 

H —> invis. 

+5% 

-0.01% 

H^qq 

+5% 

-60.05 % 

H->WW* 

+5% 

-0.18% 

H-i-ZZ* 

+5% 

-0.30% 

H —> x+x“ 

+5% 

-60.60% 

H —> yy 

+5% 

-60.79% 

H->Zy 

+5% 

-0.74% 

H —> WW* —> xvxv 

+5% 

-0.98% 


Table 6 Biases in the extracted H(Z —> qq) cross section for cases 
where the Higgs BR to a specific final state is increased by 5 %, i.e. 
BR(H X) BR(H -> X) -h 0.05. 


ified (non-SM) branching fractions and the total (visible H- 
invisible) cross section is extracted as before (assuming the 
SM Higgs branching ratios). Tab. 6 shows the resulting bi¬ 
ases in the extracted total cross section for the case when 
a BR(H ^ X) BR(H -> X) ±0.05. In all cases, the re¬ 
sulting biases in the extracted total HZ cross section are less 
than 1 %, which should be compared to the 1.8% statisti¬ 
cal uncertainty. These variations represent large deviations 
from the SM which would be observable in studies of exclu¬ 
sive final states. For example, for an integrated luminosity of 
500 fb^*, a 5 % (absolute) increase in branching ratio would 
result in an increase of 3350 HZ events in that particular 
Higgs decay topology, including an increase of 230 events 
with either Z —and Z —e+e^ decays. Such large ef¬ 
fects would be observable at a linear collider either through 
their impact on exclusive Higgs branching ratio analyses or 
they would manifest themselves as large excesses of events 
in the event samples obtained from the Z —e+e^/p+p^ 
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recoil mass analysis. It is therefore reasonable to conclude 
that unless very large BSM effects had been previously dis¬ 
covered, the hadronic recoil mass study gives an effectively 
model-independent measurement of the HZ(Z —^ qq) cross 
section. 


4 The Hadronic Recoil Mass Measurement at the ILC 

The hadronic recoil mass study presented in this paper was 
first performed in the context of the CLIC accelerator, where 
the first stage of the machine was assumed to operate at 
= 350 GeV. The study was then repeated for the ILC 
at ^/s = 350 GeV. Again a full Geant4 simulation of the 
detector response and a full reconstruction of the simulated 
events was performed. Since both studies used the same sim¬ 
ulation and reconstruction software, only small differences 
in precisions on a {HZ) from the hadronic recoil mass mea¬ 
surement at the ILC and CLIC are expected. There are two 
main effects. Firstly, because of the smaller beam spot at 
CLIC, the impact of beamstrahlung is greater than for the 
ILC, leading to a larger number of events towards lower 
values of at CLIC compared to the ILC, where is 
the effective centre-of-mass energy of the colliding electron 
and positron after the radiation of beamstrahlung photons, 
although the difference is not large at y/s = 350GeV. Sec¬ 
ondly, the ILD detector concept used for the ILC studies has 
more complete calorimeter coverage down to low polar an¬ 
gles than the CLIC_ILD detector concept used for the CLIC 
studies. Both effects will tend to degrade the hadronic recoil 
mass reconstruction for the CLIC configuration compared 
to the ILC. However, the impact is not large, as can be seen 
from Fig. 9. 




P(e ,e+) 


^^invis. 

Aa(BZ) 

CLIC 500 ft-' 

0,0 

±1.71% 

±0.56% 

±1.80% 

ILC 

500 ft-' 

0,0 

±1.57% 

±0.48% 

±1.63% 

ILC 

350ft-' 

-0.8,-f0.3 

±1.68% 

±0.52% 

±1.76% 


Table 7 Summary of the statistical precision achievable on clHZ) 
from the hadronic recoil mass analysis at = 350 GeV for CLIC and 
the ILC. The ILC numbers are shown for both zero and the nominal 
beam polarisations. 


Tab. 7 compares the statistical precision achievable at a centre- 
of-mass energy of = 350GeV for; 500fb^' at CLIC with 
unpolarised beams; 500 fb^* at the ILC with unpolarised 
beams; and 350 fb^* at the ILC with the nominal ILC beam 
polarisations ofP(e^,e+) = (—0.8,-1-0.3). For the same in¬ 
tegrated luminosity and unpolarised beams, the precision 
achievable at the ILC is approximately 8 % better than that 
at CLIC, reflecting the slightly better recoil mass resolution 






m,JGeV 


Fig. 9 The reconstructed hadronic recoil mass distributions for events 
passing the preselection cuts a) for the clear two-jet topology of the 
invisible Higgs decay analysis and b) for the visible Higgs decay anal¬ 
ysis. In both cases the distributions compare the CLIC and ILC simu¬ 
lations for 500fb^* at y'i = 350GeV, with unpolarised beams. 

at the ILC seen in Fig. 9. Since the instantaneous luminosity 
at the ILC is expected to scale with the Lorentz boost of the 
colliding beams 7e, the time taken to accumulate 350 fb^* of 
data at ^/s is comparable to the time required for 250 fb^ * at 
= 250GeV. Hence, for the nominal ILC beam polarisa¬ 
tion off’(e^,e+) = (—0.8,-1-0.3), the statistical precision of 
1.8% achievable on the HZ cross section at = 350GeV 
using Z —^ qq is directly comparable to the statistical pre¬ 
cision of 2.6% [5, 8, 9] achievable with 250fb^' of data 
at y/s = 250GeV using Z —> decays. This conclusion 

weakens the motivation for operating a future linear collider 
significantly below the top-pair production threshold. 


5 Centre-of-mass Energy Dependence of the Hadronic 
Recoil Mass Analysis 

The hadronic recoil mass analysis described above for y/s = 
350 GeV was repeated for CLIC at yG = 250 GeV and yG = 
420GeV. In each case a full set of SM model background 
processes was generated using the Geant 4 simulation of 
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Fig. 10 The two-dimensional distributions of the reconstructed Z mass, mqq, versus the reconstructed hadronic recoil mass, mrec, in the visible 
Higgs decay analysis, broken down into HZ(Z —> qq) signal and the SM background for CLIC operating at = 250GeV, y'J = 350GeV and 
,/s = 420 GeV. All events passing the preselection cuts are included. 


the CLIC_ILD detector concept. Because the complete sim¬ 
ulation of the CLIC beam is not available for these centre-of- 
mass energies; the 250 GeV samples used that same j^/s 
distribution as for i/s = 350 GeV, whereas the 420 GeV used 
the \[s' j^s for the 500 GeV option for CLIC. The analysis 
described in Sect. 3 was repeated at each centre-of-mass en¬ 
ergy using the appropriate distributions for the likelihood 
function. The binning and range used for mjec in the two- 
dimensional distribution of mqq versus mrec was optimized 
for each centre-of-mass energy. The resulting sensitivities 
are listed in Tab. 8. Compared to ,/s = 350GeV, the over¬ 
all sensitivity for A (7 (HZ) is worse at both = 250 GeV 
and = 420 GeV, although for two different reasons (ex¬ 
plained below). 

Fig. lO shows two-dimensional distributions of mqq versus 
mrec, broken down into signal and background for the three 
centre-of-mass energies considered. Sect. 3.3. For all centre- 
of-mass energies, the most signihcant backgrounds are from 
e+e^ —^ qqqq and e+e^ ^ qqfv. The qqfv background (pre¬ 
dominately from e+e^ —> WW) accounts for the broad band 
of events on the left-hand side of the background plots. This 
event population is well separated from the signal region. 
The more signihcant background arises from the qqqq hnal 
state, populating the regions with mqq ~ mz. In this region, 
the qqqq background arises primarily from e+e^ —WW, 
e+e^ —ZZ* and e+e^ —> Zy*, where the indicates an 


off-mass-shell particle; the component from e+e ^ ZZ, 
where both Z bosons are on-shell is largely suppressed by 
the preselection cuts. The board recoil mass distribution for 
the preselected qqqq background is pushed towards the kine¬ 
matic limit due to two main effects: i) the pair of jets with 
the invariant mass closest to mz is used to calculate the 
four-momentum of the assumed Z boson, in the case of the 
e+e^ —> WW qqqq background, this can lead to pairing 
of two jets from different W-boson decays; ii) for events 
with signihcant ISR or beamstrahlung, the calculated recoil 
mass (which uses the assumed centre-of-mass energy y(s, 
rather than sTs') is higher than the invariant mass of the re¬ 
coiling system. 

From Fig. 10 it can be clearly seen that the width of the 
recoil mass distribution for HZ(Z ^ qq) events increases 
with increasing centre-of-mass energy. This can be under¬ 
stood from the expression for the recoil mass: 

= (v^-Fz)^- (-Pz)^ 

= s-2Vs£'z+£'z-Pz 

« s-|-m| -2y(s(£'i -f £ 2 ), 

where E\ and £2 are the energies of the two jets forming the 
reconstructed Z boson and assuming £| — p| » m|, which 
is true for the signal region. Propagating the errors on the jet 
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Machine 



(7(HZ) 

^^vis. 

^ ^invis. 

Zl (t(HZ) 

CLIC 

250 GeV 

500 ft-' 

I36ft 

±3.63% 

±0.45 % 

±3.65% 

CLIC 

350 GeV 

500 ft-' 

93 ft 

±1.71% 

±0.56% 

±1.80% 

CLIC 

420 GeV 

500 ft-' 

68 ft 

±2.42% 

±1.02% 

±2.63% 


Table 8 Summary of the statistical precision achievable on ^(HZ) from the hadronic recoil mass analysis at CLIC for = 250 GeV, = 

350GeV and y'J = 420GeV. In each case unpolarised beams were assumed. 


energy measurements, (7i and 02 , implies that 






of+ 0 , 


Therefore, the recoil mass resolution is expected to worsen 
with increasing centre-of-mass energy due to both the yfs 
dependence and the fact that the absolute uncertainty on 
the jet energies increases with jet energy (Og ^ O.OSis) and 
therefore with centre-of-mass energy. The increasing width 
of the recoil mass distribution accounts for the increase of 
4<Tinvi.s. with ^/s, listed in Tab. 8, and the larger value of 
ZiOvis. at y/s = 420GeV. However, despite the better recoil 
mass resolution, the sensitivity to ZiOvis. at ^/s = 250 GeV 
is significantly worse than for the other centre-of-mass en¬ 
ergies considered. The reason for this can be seen clearly 
in Fig. 10. At yfs = 250GeV, HZ production is not very 
far above threshold and the recoil mass distribution is rela¬ 
tively close to the kinematic limit. This is the region popu¬ 
lated by the large qqqq background passing the preselection 
cuts, resulting in a greatly reduced separation between sig¬ 
nal and background in the variable that provides the best 
distinguishing power, namely ■ 


6 Summary and Conclusions 

This paper presents the first detailed study of the potential of 
the hadronic recoil mass analysis at a future linear collider, 
both for visible Higgs decay modes and possible BSM in¬ 
visible decay modes. By combining the analyses for visible 
and invisible modes, it is shown that the measured e+e^ 
HZ(Z —?> qq) cross section does not depend strongly on the 
nature of the Higgs boson decay and thus provides a model- 
independent determination of gnzz- The statistical precision 
achievable at CLIC operating at ^/s = 350 GeV with 500 fb^ * 
of data with unpolarised beams is zIcThz « ±1.8%. A sim¬ 
ilar precision is obtained for the ILC with 350fb^^ and the 
nominal beam polarisation ofP(e+,e^) = (-1-30%,—80%). 
In both cases the branching ratio to invisible decay modes 
can be constrained to BR(H —invis.) < 1 % at 90 % confi¬ 
dence level. It is demonstrated that yfs = 350 GeV is likely 
to be close to the optimal energy for the hadronic recoil mass 
analysis; at lower centre-of-mass energies there is less dis¬ 
crimination between signal and background and at higher 


centre-of-mass energies the measurement is limited by the 
worsening recoil mass resolution. 

It is often stated that operation of a future e+e^ linear col¬ 
lider close to threshold {^/s ^ 250GeV) is necessary to pro¬ 
vide an absolute measurement of the coupling between the 
Higgs boson and the Z boson, gnzz- This is based on the de¬ 
termination of gHZZ from the recoil mass analysis for e+e^ —>^ 
HZ(Z —The results presented in this paper show 
that, for a comparable running time, a statistically more pre¬ 
cise measurement can be obtained from e+e^ —>■ HZ(Z —>■ 
qq) events at ^/s = 350GeV. This conclusion argues against 
initial operation of a future linear collider at significantly be¬ 
low the top-parr production threshold. 
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